arXiv:1503.04820vl [hep-th] 16 Mar 2015 


NORDITA-2015-29 


Baryon and chiral symmetry breaking in holographic QCD 

Alexander Gorsky/’^ Sven Bjarke Gndnason,^’"^ and Alexander Krikun"^’® 

^Institute of Information Transmission Problems, 

B. Karetnyi 19, Moscow, Russia 
^Moscow Institute of Physics and Technology, 

Dolgoprudny IfllOO, Russia 

^Institute of Modern Physics, Chinese Academy of Sciences, 

NanChangLu 509, Lanzhou 730000, China 
^Nordita, KTH Royal Institute of Technology and Stockholm University 
Roslagstullsbacken 23, SE-106 91 Stockholm, Sweden. 

^Institute for Theoretical and Experimental Physics (ITEP), 

B. Cheryomushkinskaya 25, 117218 Moscow, Russia 

Abstract 

We study the relationship between chiral symmetry breaking and baryons in holographic QCD. We 
construct a soliton with unit baryon charge in the presence of a nonzero mean value of the scalar 
bifundamental field, which is dual to the chiral condensate. We obtain a relation between the chiral 
condensate and the mass of the baryon and find in a clear-cut way that at large values of the condensate 
the holographic soliton is no longer located on the IR wall. Instead it is split into two halves, which 
are symmetrically located on the left and right flavor branes. On the other hand we find that the local 
value of the quark condensate is suppressed in the core of the soliton, which is evidence for a partial 
chiral symmetry restoration inside the baryon. 
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I. INTRODUCTION 


Baryon physics in holographic setups has been addressed since the time the hrst holographic 
QCD models have been developed. Th^r^is a wealth of works related to the interpre^Bon of 
the baryon as a holographic instanton 


even high-density baryonic matter 


-yl , the QCD currents sourced by the baryon [4-l6| and 
7|-l9|]. Unfortunately though, the studies of the interplay 
between baryon physics and physics of chiral symmetry breaking is undeservedly rare: Ref. 101 
is the only example that we are aware of. But from the other hand these relationships remain 
one of the longstanding questions in QCD |n| and should be dehnitely addressed with the new 
tools provided by the holographic duality, which proved to be useful in describing both chiral 
symmetry breaking and baryons in QCD. 


Usually the AdS/QCD model, whic 
of) the Sakai-Sugimoto model (SS) 


h is used for the studies of baryons, is (a modihcation 
13| . This model is based on the system of Nf D8 — D8 
branes, which are embedded in the curved hyperbolic space produced by a stack of Nc DA 
branes. One of the dimensions, r, is compact and serves for breaking supersymmetry. The 


compactihcation radius of r vanishes at some point and thus the background space has the 
shape of a cigar |l^. The location of the tip of the cigar plays the role of the Aqcd scale. From 
the holographic point of view, the tip serves as a boundary beyond which nothing can propagate 
as the space practically ends there. The D8 — D8 branes embedded in this background are 
located at the opposite (or almost opposite |l^) points of the r-circle and consequently they 
merge together when the circle collapses. The gauge held. A, which is a massless mode of the 
open string with two ends connected to the D8 branes, belongs to the adjoint representation of 
the havor group U{Nf) and is dual to the quark current operator of QCD. One can introduce 
the coordinate U along the branes with the origin on the tip of the cigar, which goes from — oo 
on the D8 brane to -|-cxd on the D8 brane. The values of the gauge held A at —oo may be 
thought of as being coupled to the left quark current = q^-^'yA°‘q while the values at oo 
are coupled to the right current = q^-^'^A'^q- Accordingly one can call the D8 branes 
Left and the D8 branes Right. 


The baryon in the Sakai-Sugimoto model is an instanton of the held A [l|. The baryon 
charge equals the topological charge of the instanton because the baryon current is sourced in 
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the Chern-Simons (CS) term of the action by a topological charge density. One can hnd that 
the instanton in the SS model is well approximated by the usual flat BPST instanton. This 
holds because the size of the instanton is parametrically small; the center of the solution is 
located on the tip of the cigar and one can neglect the curvature of the metric at the scale of 
the solution radius. The vanishing difference between spatial (xj, i = 1) • • •) 3) and holographic 
(U) parts of the metric leads to an approximate S'0(4) symmetry group of the space [^. In 
this setting it is very convenient to use the usual radial ansatz of the BPST instanton with the 
radial coordinate = U'^ + The solution is then located at p = 0. 

For our purposes, though, the Sakai-Sugimoto model is not entirely satisfactory. The obstacle 
is that it is not at all easy to describe chiral symmetry breaking and the associated chiral 
condensate in this setup. This problem was analyzed in Refs, la, llTI] and it was shown that 
the holographic held dual to the scalar quark current qq is a tachyonic mode of the string 
stretched between the D8 and the D8 branes. The tachyonic held condenses on the tip of the 
cigar, thus providing a nonzero vacuum expectation value (VEV) of qq, the chiral condensate. 
There are, however, several complications in this setup. The construction of the tachyonic 
DBI action is ambiguous, the backreaction of the condensed open string on the geometry of 
space should be taken into account and so on. On top of that, when constructing an ehective 
5-dimensional holographic model for mesons it is impossible to treat the open string as a local 
held. 

The chiral symmetry breaking is relatively easy to implement in the other, “bottom-up”, 
approach of holographic QCD. Instead of constructing the brane system which would reproduce 


the ehective theory of mesons, one builds the mode’ 


operators, in the 5-dimensional AdS background 18 


a including the helds, dual to the QCD 
. In the “hard-wall” (HW) model 18 1 


the space ends on the wall located at a hnite value of the holographic coordinate Zm- This wall 
(we will call it the IR wall) serves the same purpose as the tip of the cigar in the SS model 
by breaking conformal symmetry and providing the scale Aqqd- There are two gauge helds L 
and R, which are dual to left and right quark currents. The held, which is dual to the scalar 
quark operator qq, is just a bifundamental scalar X. Its vacuum prohle is described by the two 
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parameters 


1 . ox 

^0 = + az ), 


( 1 ) 


which W the holographic dictionary 20j are related to the quark mass m and condensate (see 


Refs. 


21 


231] and the derivation in Sec. IVII) 


(qq) = 


27r^ 


a. 


( 2 ) 


In the model there is no internal mechanism, which would £x the value of the chiral condensate. 
Instead one obtains a nonzero value of a by imposing a Dirichlet boundary condition on X at 
the IR wall 


X 


1 / S N 

= -(mZm + (Tz^j = const. 


( 3 ) 


One can obtain the correlation functions of various operators in the setup with nonzero prohle 
of the bifundamental scalar and study the dependence of various observables on the chiral 
condensate, which is just a parameter of the model. Unfortunately, the treatment of the baryon 
in this setup is not as straightforward as it is in the SS model. There are two gauge helds, 
which should realize the topologically nontrivial solution, and on top of that, there is a scalar 
with a nonzero VEV, which interacts with both of them. One should also pay special attention 
to the boundary conditions at the IR boundary, as they are arbitrary in the construction and 
should be additionally fixed in the model. Concerning the baryon, the question of IR boundary 
condition becomes substantially important because the corresponding topological soliton falls 
on the IR wall, a phenomenon similar to the localization of the SS-instanton on the tip of the 
cigar. The baryon physics in HW model without a chiral condensate was addressed in Ref. 24l |. 


The bifundamental scalar was included in Ref. [10|, however the dependence of the baryon 


mass on the condensate 
QCD sum rules in Ref. 
of the hard wall model 


las not been investigated. Motivated by the formula derived using the 


25 


m two of US have analyzed the baryon-mass origin in the framework 


26| and found that at large values of the chiral condensate, there 


are clear indications that it dominates. However the accuracy of the numerical calculations was 
not high enough to make precise claims. In this paper we perform a detailed analysis of the 
different aspects of the impact of the chiral symmetry breaking on the baryon solution. We 


4 

















will clearly see the internal structure of the solution and two regimes in the dependence of the 
baryon mass on the chiral condensate. 

The paper is organized as follows. In Sec. [TTl we introduce our setup, which we try to keep 
sufficiently general and applicable to a wide range of holographic QCD models. We study the 
qualitative structure of the baryon solution in the presence of the chiral scalar held in Sec. IIIII 
The quantitative support for this treatment is given by the numerical solution in Sec. IIVI 
Secs. El and El] are devoted to the calculation of the baryon mass and local mean values of the 
QCD currents in presence of the baryon. We conclude in Sec. I VI II The analysis of the precision 
of our numerics is outlined in the Appendix. 

II. THE MODEL 

It is actually possible to make a qualitative connection between the seemingly different 
setups of the Sakai-Sugimoto and the “hard wall” models. In order to build up our intuition 
we use the following “folding trick”: Consider the cigar-bended D8 brane of the SS model. 
There is a single gauge held A on this brane, but let us call it L when it is considered on the 

Left brane and R on the Right brane (see Fig.[T]). On the tip of the cigar L = R, oi course. 

Next, let us assume trivial dynamics along the r direction, which will allow us to neglect the 
separation of the branes along r altogether, ehectively “folding” the cigar-shaped brane into 
a single sheet. More specihcally, we cut the brane along the tip or f/ = 0 line. On one side 
we have the L held and U G [0; -|-oo) and on the other side the R held and U G (— 00 ; 0]. In 
order to get consistent coordinates, we perform a.U ^ —U transformation on the Right brane, 
which also means Ru —)■ —Ru- It is convenient to introduce the inverse holographic coordinate 
z = 1/{U + 2 ;“^). Then the tip of the cigar (f/ = 0) corresponds to the maximum z = Zm and 
the boundary values are assigned at z = 0. At the end of the day we have two gauge helds 
propagating on a strip 2 ; G (0,^^], which is qualitatively the same as in the HW model. The 
important additional information, that we get here, is the boundary conditions on the gauge 
helds at the IR boundary. Because we know that L and R are the same held there, we get 

IR : Lz -|- Rz = 0, Lf^ — R^ = 0. (4) 
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In this picture the tachyonic string, which was responsible for the chiral condensate in the SS 
model, collapses to the local bifundamental tachyonic scalar X. The SS baryon, which was 
originally located at 17 = 0, is now split into two halves, one composed of L and the other of 
R, which are folded together. Thus we expect that in the HW model, the baryon would appear 
as one half of the instanton lying on the IR wall. 



FIG. 1; Scheme of the “folding trick” explained in the text. 

As we see, in any case the holographic QCD model with two quark flavors Nf = 2 can 
be described by dynamics of left (T) and right {R) gauge fields of the group U{2)l x U{2)ji 
and a bifundamental tachyonic scalar X in 5-dimensional curved space. We use Hermitian 
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gauge-group generators and treat the gauge connection as Hermitian matrices^: 

L = L — L = L ^“ 1 “ , a = 1 ,..., 3 , ( 5 ) 

= y, tr(t“t^) = [t“, t’’] = 


where cr“ are the Pauli matrices. The background space can be described by the mostly-minus 
metric of the general form^ 

ds^ = h{zYdt^ — h{zYdxl — k^z^dz"^, t = xq, z = i = 1,... ,3. ( 6 ) 

For instance, in the “hard-wall” model k{z) = h{z) = Ljz] in Sakai-Sugimoto model 121 
h{z) = {L/zf/\ k{z) = {L/zy/\l - L is a curvature scale and the holographic 

coordinate 2 ; is bounded from above: 2 G (0; Zm]- The Yang-Mills action is 

Sym = I d^xdtdzi^ - -^TT^k{z)FLl - 2k{z)FLl + 2 ^ - F^l,) + {L ^ R)) ( 7 ) 

h{z'^^ 


+g\Tilh{zf (FoX^-AX^) 


k{z] 


-{D^Xy - h{zfk{z)mi^\X\ 


where the covariant derivative is = d^X — iL^,X + iXR^ and the g auge coupling 


23l | or a specihc 


and normalization of the scalar gx are to be hxed by phenomenology 0,121^ 
top-down construction 0, 16]^. The mass of the scalar mx is dehned by the dimension of the 
corresponding operator and equals rrix = —j 2 - 


The Chern-Simons term is present as well 


Scs = 


A 

247r2 


- <’LATr(AAA)-(Lo A) 


( 8 ) 


One can immediately spot that the non-Abelian part of the gauge held, which has nonzero 
topological charge in the 4-dimensional spacelike slice, provides a source for the temporal com¬ 


ponent of the Abelian held, which is dual to 


baryon charge of a given held conhguration 


0 


he baryon current. Thus we can identify the 


24| 


Qb — 


1 




d^xdz i^pamnpalk 


^amn ^alk\ 

R R J "> 


R 


m = 1 ,..., 3 , 5 . 


( 9 ) 


^ The corresponding definition of the gauge field strength is F^i, = dfiA^ — — i[A^,Ai,] 

2 Note the peculiar numbering of the coordinates: we never use Xi as it usually stands for Euclidean time and 
our treatment is Minkowskian. For the holographic coordinate we use X 5 instead to avoid possible confusion. 
^ Importantly, in the Sakai-Sugimoto model, is inversely proportional to the ’t Hooft coupling of the dual 
field theory, which is assumed to be large. 
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In order to construct the solution in the 4D spacelike {xi, 2 ;)-plane we adopt the cylindrical 


ansatz proposed in Refs. 


i],E, 


271], taking r = xf and 2 ; as the cylindrical coordinates. The 


4D gauge potentials are parametrized by all possible tensor structures which link the spatial 
and SU{2) group indices. It is convenient to work with the vector and axial combinations, 
dehned as 

L = V + A, R = V -A. (10) 

In order to preserve 3D parity, we choose the P-odd and P-even tensors for the spatial com¬ 
ponents of Vi = {Li + Ri)/2 and Ai = {Li — Ri)/2, respectively. Similarly, the time and 2 ; 
components should be parity even for V and odd for A. Hence we obtain the following ansatz 
for the gauge helds 


T/a_ l + V 2 {r,z)^ Xk 


^ 5 “ = 0, 




XjXa 


For the scalar we have 


+ A{r,z) 


XjXa 


Xn 


At = Az{r, z) — , 


Vo = v{r,z), 

Aq = 0 . 


( 11 ) 


X = xi(r, z)- + IX 2 V, z) -. 

2 r 

The action ([7]) with this ansatz is S = f dtdrdzi^CyM + Rui + '^cs)- 
271 


Rym — — 


9l 


- drA^]‘ 


271 g\{ r 


k{z) 

2 

k{z) 

Mzf 


k{zfVrril h{zf'D^ril 


k{z 

g HZA 

k{z) 


k{zfVrxl + h{zfv,xl + m\h{zfk{zfxl 


[X2V2 + Xlhl]' + (1 - dlf 

9x95 ' 


a = 1,2, 


Cui — I 


dvr ^ '"2 r 


2^i/ k{z) 


h{zY{dzvY + k{zY{drVy 


k^cs = 


El 

277 


2e°‘^V^r]o,Vrr]^ + {d,Ar - dr A,) {l - vl) 


( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


where the Abelian covariant derivative is = 9*?7q, -|- A^,e°‘^ri^, and similarly for Xa- We see 
that the problem boils down to the 2D Abelian Higgs model with two complex scalars rj and 
y, which interact with each other, plus some potentials as well as interactions with the real 

























scalar v, which has the opposite sign of the kinetic term. One can expect to have vortices as 
topologically nontrivial solutions in this model 271. 


III. STRUCTURE OF THE SOLUTION 


In order to study the structure of the soliton it is useful to introduce the phases of the 
complex scalars 


+ »')2 = rje'". Xi + iX2 = Xe’’’- 


(19) 


First we note that the baryon charge ([9]) can be rewritten in the form 


Qb — — 

TT 


drdz 


dr{'Dxricrife“^) - + (diA - drA,) 


( 20 ) 


Once we demand that, in a particular gauge, the first two total derivative terms vanish^, we 
are left with the standard topological charge of Abelian Higgs model. Assuming the pure gauge 
condition Aj = did + const at the boundaries, we realize that the topological charge equals the 
winding of the phase 6 around the boundary of the spatial patch. 


Qb — 


1 

TT 



■ dr6 


drdzB). 


( 21 ) 


One should note though that the baryon, which is a solution with Qb = 1, should have the 
phase 9 winding by vr. Hence the baryon in our model is a half-vortex.®. The only way a 
half-vortex can be realized as a smooth solution is by having its core located on one of the 
boundaries. There is only one boundary of the [0, oo)^ x {0,Zm\z patch which is suitable for 
this. Indeed, in the core of the vortex the modulus rj should necessarily vanish. The second 
potential term in Eq. (1T6|) does not allow rj to deviate from 1 either at r = 0 or at z = 0 where 
k{z) is singular. The core cannot be located at infinity as well, because this will require a 
nonzero derivative drfj at r —)■ cxo and lead to a divergence in energy due to the rj kinetic term 


^ The boundary terms at r —>■ 0, r —>■ oo and 2 —>■ 0 vanish because the finite energy condition forces 'Dz^r'Ha = 0; 

while the term at z = Zm vanishes because of the boundary condition r]i\^ =0. See details below. 

^ This result is consistent with our general intuition of “folding” the Sakai-Sugimoto model, where the baryon 
is a full vortex, which lies on the tip of the cigar. 
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ffT^ . Thus we see that the core of the rj half-vortex should be located on the IR wall 2 ; = Zm, 
which corresponds to the tip of the cigar in the Sakai-Sugimoto picture. 

At the next step we should consider the interaction of the scalar helds rj and x- h is governed 
by the hrst term in Eq. flTbl) . which conveniently can be represented in terms of the phases flT^ : 

[X2V2 + XiVi ]' = cos(0 - 7)2. (22) 

We see that this interaction effectively “locks” the phases 6 and 7 , requiring 6 — j = ^ + irn. 
Therefore, if 6 has a winding vr, the same holds true for 7 and the scalar x realizes a half¬ 
vortex as well. The location of the core of this vortex is different though. At inhnity we need 
to connect the baryon solution with the vacuum, where the prohle of the X held is given by 
the quark mass and condensate ([T]), thus the modulus x cannot vanish. The core cannot be 
located at 2 ; = 0 as well, because according the holographic dictionary the boundary value of 
X is dehned by the source of the corresponding operator, which is the quark mass. On the IR 
wall, X is fixed by the boundary condition (|3]) which defines the quark condensate®. At the 
end of the day, the only remaining place for the core of the x half-vortex is the center of the 
soliton r = 0. 


® We note 
behavior 


here, that this condition is absent in the Sakai-Sugimoto model, 
of the tachyonic open string near the tip of the cigar 16| . 


Instead one should consider the 
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FIG. 2: Schematic structure of the soliton. Vr;,V;^ - half-vortices of i] and x scalars. Arrows show 
the effective forces produced by the metric {h{z), k{z)), kinetic terms of x and {IC^,ICa) and self 
interaction with the Abelian field in the CS-term. 

At this point we see that the holographic baryon is a soliton consisting of two interacting 
half-vortices of two scalar fields. One is located on the IR wall and another at the center of the 
soliton (see Fig. [2]). Now it is important to consider the energy of the solution and ensure that 
there exists a mechanism which stabilizes its size, preventing the soliton from collapsing. 

It was pointed out in Ref. [l| that due to the metric factors h{z),k{z) the energy of the 
soliton is roughly inversely proportional to its position. Hence the solution tends to lie at the 
largest possible 2 ; = Zm- In other words, the metric provides a force, which pulls the solution 
to the IR wall. 

On the other hand, due to the IR boundary condition on X, the core of the x vortex cannot 
approach IR wall too close. This would lead to a large gradient of x, which evolves from zero 
in the core to x ~ cr on the wall. Hence the kinetic term of x and IR boundary condition 
(l3l) provide a counterforce, which repels the solution from the wall. 

One can see though that this counterforce is proportional to the radius of the solution 
squared. Therefore it becomes negligible when the size of the soliton is small. Moreover, the 
whole X energy term flTSl) is proportional to r^, so it contributes (in addition to the gauge held 
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to the pressure, which would shrink the solution to zero size. 


24| that the important counterforce to this shrinking, which stabilizes the 


kinetic term fll3J) 

It is known [1 

radius of the solution, is provided by the interaction of the soliton with the Abelian gauge held 
via the Chern-Simons term ffTSj) . The smaller the radius of the solution becomes, the higher the 
density of topological charge is, which sources the temporal component of the vector Abelian 
held, and the larger its gradient is. This results in the ehective internal pressure which hxes 
the size of the soliton (see Fig. [2]). 

This treatment shows an interesting ehect, which is produced by the chiral scalar held in the 
“hard-wall” model. Once the soliton radius is hxed by the CS term, the IR boundary condition 
(]3|) together with the x kinetic term (IT^ produce the force which repels the soliton from the IR 
wall. The larger the chiral condensate is, the stronger the repulsion is. At large enough values 
of a, this repulsion will dehne the equilibrium z coordinate of the soliton center and therefore, 
will directly ahect its energy. Hence one should expect a substantial dependence of the baryon 
mass on the chiral condensate. In the Sakai-Sugimoto model, though, the radius of the soliton 
is inversely proportional to the ’t Hooft coupling and is therefore parametrically small l|, l2|. 
This renders the repulsive force parametrically small as well and consequently it is questionable 
whether in the Sakai-Sugimoto model any impact of the chiral condensate on the baryon mass 
can be observed. 


IV. NUMERICAL SOLUTION 


In order to conhrm the expectations outlined above we construct the soliton solution nu¬ 
merically. For the numerical simulation we choose the “hard-wall” model, so the metric is 
unperturbed AdSs: h{z) = k{z) = Ljz, and the scalar asymptotic prohle is given by Eq. ([T]). 
The couplings in the model have the values 3, 3, 23] 


9l = 


127r^ 


9l 


(23) 


^ It was noted in Ref. [l^ that the solution can be stable even in the absence of the Chern-Simons term. Our 
treatment shows that this should not be the case. Indeed, when the solution shrinks the contribution from y 
weakens and cannot prevent further shrinking. 
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The IR wall is located at Zm = (323 MeV)“^ which is hxed by the mass of the p-meson 1^. The 
parameter a is related to the quark condensate as in Eq. ([2]). In what follows, unless explicitly 
stated otherwise, we will measure all the dimensional quantities in terms of Zm and rescale the 
curvature radius L = 1 . 

Because we are dealing with a gauge held theory, hrst of all we need to hx the gauge. It is 
useful to adopt the Lorenz gauge dzA^ + drA^ = 0. The gauge hxing condition is enforced by 
adding the Lagrange term to the action 


A [dzAz + drArf . (24) 

With this choice of gauge hxing, the equations for Az and Ar become elliptic and the boundary 
value problem has unique solution. The resulting physical solution should not depend on the 
value of A. 

Next we need to choose the boundary conditions consistently with the given gauge. In the rj 
half-vortex the phase 6 changes by tt along the boundary of space. We can chose the gauge in 
such a way that the total winding occurs along the r —)■ cx) boundary, where 6 changes linearly 
^lr-j.oo “ constant everywhere except the jump by vr in the core of the r] 

half-vortex on the IR boundary. This choice of 6 ensures us that 772 = —1 at 2 ; = 0, hence there 
are no source terms for the axial current, 772 = —1 at inhnity, hence the currents decay there, 
and 771 = 0 on the IR wall which is required by the IR condition (jl]). Due to the interaction 
between the scalars y and rj fl 22 |) . the phase difference {9 — 7 ) should be either — | or | on 
all the boundaries. Therefore the phase 7 winds linearly along the r ^ 00 boundary as well: 
l\x=i ~ constant along the z = 0 and z = Zm boundaries, having a vr jump in the 

core of the y half-vortex at r = 0. One can check that this choice of phase guarantees the 
constant value of the quark mass ?n at 2 ; = 0 . 

The pure gauge condition at r —)■ 00 will then hx Az to be a constant: Az\^^^ = dz9 = 
dz') = Hence according to the Lorenz gauge, we get drAr\^^^ = —dzAz = 0. As we 

discussed in the introduction, on the IR boundary Ar _ =0 due to the merging of the Left 

and Right branes at this point (|1]).® This does not apply though to the z component and thus 


In Ref. iioi it was also noted that this condition is required by gauge invariance of the CS term in the 
“hard-wall” model. 
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the gauge condition will take the form _ = —drAr = 0. The boundary condition on 

the UV boundary z = 0 is dictated by the absence of the sources for the spatial axial current, 
hence Ar|^_Q = 0 and dzA^ = —drAr = 0. Finally in the center of the soliton the pure gauge 
condition states Az\^_^ = 0 and consequently drAr\^_^ = 0. The boundary conditions for the 
Abelian vector field v are set by the absence of sources, regularity at the center and on the IR 
wall as well as the fall-off condition at infinity. Reexpressing these results in terms of the scalar 
fields rja, Xa we get the set of the boundary conditions shown in Table [H which we will use for 
the numerical solution. Conveniently, all of them are either of Dirichlet or of Neumann type. 


r —)■ cx) 

z = 0 

r = 0 

Z = Zm 

^i=cos(-f+vr^) 

r?i =0 

r/i =0 

r?i =0 

r? 2 =sin(-f+7r^) 

??2 = -l 

??2 = -l 

dzV2 = 0 

Xi = (mz + az^) cos (vr A-) 

Xi = (mz + az^) 

drXi = 0 

Xi = — (mz + az^) 

X 2 = (mz + az^) sin (vr^) 

'■ ' \ Zm' 

X2 = 0 

X2 = 0 

X2 = 0 

o 

II 

Ar = 0 

o 

II 

cd 

Ar = 0 

= zJt: 

O 

II 

CO 

A, = 0 

O 

II 

CO 

f = 0 

v = 0 

O 

II 

o 

II 


TABLE I: Boundary conditions for the fields used in the numerical calculation. 


We look for a solution to the equations of motion, which follow from Eqs. (lT3|) . (ll71fT8i) . In 
order to reduce the problem to the calculation on a square patch [0, l]^; x [0, l]y we introduce 
the rescaled coordinates 


r = c tan 



Z = ZmV, 


(25) 


where the dimensional constant c specifies the radial scale on which the solution will be best 
resolved. On this patch the homogeneous grid x Ny is introduced and we construct 7 
equations (one for each dynamical held) on each node by the nearest neighbor discretization of 
the derivatives in the equations of motion. 

Special attention should be paid to the boundaries r = 0, z = 0 and r —)■ oo, as the equations 
of motion are divergent there. We do not perform the regularization by stepping out from the 
boundary. Instead, for the helds with Neumann boundary conditions, we expand the equations 
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of motion near the singular point and solve the leading order contribution. This allows us to 
impose the boundary conditions directly on the boundaries. 

After discretization we obtain 7 x {N^ — 2) x {Ny — 2) algebraic equations in internal points 
of the grid plus + 4:Ny equations on the boundaries and solve the resulting system both 
by Newton-Raphson method®and the relaxation method, independently, in order to have 
a good check of the numeric results. The analysis of numerical accuracy is described in the 
Appendix. 

As a starting held conhguration we take the superposition of two vortices according to Fig. O 
which has unit topological charge ([9]). At the end of the day we were able to hnd solutions 
with a number of different parameters a and m on the grids 60 x 60, 120 x 120, 240 x 240 and 
observe good convergence and matching between the two methods (several other grids were 
used as well, see the Appendix). The procedure is run until the change in the held values on 
each step falls below 10“®. For the solutions we check that they do not depend on the value 
of the Lagrange multiplier A, ensuring that the gauge hxing condition fl2Tll is fulhlled. On top 
of that, we check that the observables do not depend on the choice of the parameter c, which 
dehnes the radial coordinate transformation fl25|) . 


V. BARYON MASS 


Before we proceed with the analysis of various baryon parameters, it is interesting to look 
at the general behavior of the solution. The distribution of the topological charge at various 
a is shown in Fig. |3l At very low a the solution is concentrated near the IR boundary and 
does not feel the presence of the scalar held. At intermediate values of the chiral condensate 
the charge distribution is pulled downwards by the scalar “repulsion force” described above. 
Finally, at large a one can see that the solution is detached completely from the IR wall and 
is localized at intermediate values of the holographic coordinate 2 ;. Another observation is the 

® Strictly speaking we use the two-step procedure described in Ref. [^: first we solve for all the fields except v, 
then we solve the linear equation for v and so on. This is due to the inverse sign of the v kinetic term, which 
prevents us from finding the solutkm as a simple minimum of the action functional. 

We use Wolfram Mathematica 9 [28| for deriving the equations and compiling the numerical code and the 
LinearSolve [] procedure therein to invert the resulting matrix. 
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dependence of the radial coordinate of the solution on the value of the condensate. At large 
(j the shrinking force produced by the scalar part of the action grows and it is harder for the 
CS term to stabilize the radius of the solution. All these observations are in perfect agreement 
with the intuition developed in Sec. IIIII 


(t‘^^= 0.22 1.35 cr‘''3=2.93 



0 0.97 oo 0 0.84 oo 0 0.41 

FIG. 3: Topological charge density of the soliton at various values of the scalar field boundary value. 
At large a the soliton is detached from the IR wall. The radial (horizontal) coordinate and a are 
measured in units of Zm = (323 MeV)“^. 


In order to study the mass of the holographic baryon we hrst regularize the energy by 
subtracting the inhnite contribution from the vacuum conhguration of the scalar held ([I]): 

2r^ 


^reg 


dzd^x 


—C + 27rgx—^ {rn? — 


The mass of the holographic baryon for different values of a is shown in Fig.H] (the value of 
quark mass here is m = 0.01z“^ = 3.23 MeV and we use Nc = 3 in order to obtain the numeric 
results). As anticipated, at low values of a, when the solution lies on the IR wall, the mass is 
governed by the conhnement scale. In the opposite case of large a, the position of the solution 
is only controlled by the scalar “repulsion force” and hence its mass is directly proportional to 
the chiral symmetry-breaking energy scale The intermediate region, where the effects 

of chiral and conhnement scales are comparable, is quite narrow, so using Eq. (|2]) we may 
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approximate the result for the mass of the baryon as 


Mb = NcZ^ max 


1/3 


= max 


0.92;-0.42+ 0.97 
887;-407+ 940 MeV. 


It is important to check that the N^. scaling of the baryon mass is linear due to (qq) N,. 
Interestingly, the physical value of the baryon mass Mb = 940 MeV corresponds to the point 
where the scalar repulsion just becomes signihcant (as claimed in 26|). 


^BZm 



FIG. 4: Dependence of the baryon mass on the chiral condensate. The solid lines show the asymptotes 
at small and large condensate values. The dashed line shows the physical value of the proton mass. 


VI. MEAN VALUES OF QUARK OPERATORS 


Now we calculate the mean values of the currents induced by the baryon. The vacuum 
expectation value of the operator is obtained by taking the variation of the action on the 


classical solution with respect to the boundary value of the corresponding held 


20 


29j. Given 


the held content of the model, we are able to calculate the mean values of vector and axial 
quark currents as well as the scalar quark bilinear. It is important to note also that the proper 
formulation of the operator/held duality is possible only in a specihc gauge in the 5D theory, 
namely in the axial gauge = 0. Our numerical solution was obtained in diherent, Lorentz 
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gauge (1211), hence a gauge transformation should be performed before we can obtain the mean 
values. Let us denote the helds in the axial gauge with bars: = 0. The action flTdll evaluated 

on the classical solution in axial gauge can be recast as a boundary term^^ 

Sci = j dtdr 47 rr^ (~^) ^ ^ + X 2 dzX 2 ] 

According to the equations of motion, the solutions near the boundary are represented by the 
Frobenius series 



A, = ^ ^ 


Va = + ha + • • • , 


Xa = X^a^Z + + 


(26) 


Substituting these into the action we hnd a divergent contribution as e —?■ 0 which comes 
from the x term: Sdiv ~ e“^(xa^)^. This must be regularized according to the holographic 
regularization procedure 


29l | with a boundary counter term 


Set- = 


dtdr 4?rr^ fXf L [j.i(e)^ + X2(<X] ■ 


(27) 


In order to get the mean value of the 3D spatial current, one needs to take the variation of the 
regularized action Sreg = Sd + Set. with respect to the boundary value of the corresponding 
spatial component of the 3D vector held A® or It is convenient to express these boundary 
values in terms of the 2D helds of the ansatz flTTl) . This will reduce the problem to the variations 
with respect to ria \ X« ^ and Therefore in terms of the Frobenius coefheients (126|) we get^^ 

P) 

(28) 


= 

2 


9i 

{J^)t = 

2 

I 2 


9i 

{m) PI = 

6 


^iak , 

r 


( 2 ) 

Vi 


XiXa 


T a 


( 2 ) 


Xj.Xn 


Xi 


(3) 1 , . (3) 


+ % 


P7 


2 ■ r 

The values of the coefheients can be calculated after we transform the numerical solution to 
the axial gauge. With the gauge function A, the helds in Eq. flT3l) transform as 

A = A + d\, f/o = ? 7 a cosA — e“^?7/3sin A, Xa = Xa cos A — e“^X/3 sin A. (29) 


The contribution from z = is canceled by an appropriate boundary term on the IR wall and the contribution 
from the CS term vanishes because of the zero boundary value of the Abelian vector field v. 

Note that if one substitutes the vacuum solution ([T]) into the expression for {qq) one gets the relation ([2)) 
between cr and {qq). 
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In order to reach the axial gauge Az = 0 while keeping the value of the quark mass m 


constant, we choose 


A = - 


dzA, 


(30) 


In the vicinity of the boundary, it is expanded as A(e) = —eAz{e). At the end of the day the 
coefficients in Eq. (l28j) are expressed in terms of the components of our numerical solution 




(3) _ _ a ^ 

2z " z 


= —d fi 

'la o '-'z'h 


2z 


,( 2 ) ^ 


= - dzAr 

2z 


1 

■ 2 
1 
2 
1 
2 


-Vzi — ]+e‘^^AzVzi^ 


X/3 


(31) 


'I^zVa + d^^Az VzTjp 
{dzAr - drAz) 


Given these expressions we can calculate the mean values of the currents 

It is interesting to study the local scalar VEV. In a hedgehog Skyrmion it behaves as 

{qq){r) = {qq)o ( cos(6') +i -sin(6') 


(32) 

where {qq)o is the value of the chiral condensate at inhnity, and 6{r) is a chiral phase which 
goes from tt in the center of the Skyrmion to 0 at inhnity. In the analysis of our solution we 
hnd that the phase of the scalar VEV behaves exactly as expected (see Fig. [5]), but on top 
of that, the modulus of the chiral condensate is suppressed in the center of the baryon. This 
observation agrees with the lattice studies 30| of the partial chiral symmetry restoration inside 
baryons. Spectacularly, while the value of the chiral condensate affects the baryon mass, the 
presence of thee baryon has its effect on the scalar quark density as well. 


The other interesting observable is the nucleon axial charge qa- In Refs. 


10 


3l| the value 


of the axial charge calculated for a holographic baryon was almost 2 times smaller than the 
phenomenological value qa ~ 1.25, although the other baryon observables exhibited good agree¬ 
ment. One should expect, though, the dependence of qa on the pion mass, which is related to 
the chiral condensate by the Gell-Mann-Oakes-Renner relation. Thus it is interesting to check 
whether this value changes signihcantly for different values of a. For this purpose we reproduce 


the calculation of Ref. 


10| and study the dependence of qa on the chiral condensate. If one 


parametrize the mean value of the QGD axial quark current by the following radial functions 


(J^)“(r) = A(r)V“ + A2(r): 


(33) 


19 


















{q q) 

(qqh ^ 



FIG. 5: Relative local mean value of the scalar operator (qq) and its phase 9 for various solutions 
with = {0.31, 0.37, 0.62, 1.50, 2.10, 3.10} counting from top to bottom. 


the axial form factor is expressed as 321 


GA{q^) = 


Mb\ Stt 

/ Y 


r‘^jo{qr)Ai{r) 


ji{qr)A2{r) 


and the axial vector charge qa is 

OO 

Sa = GA(q^ = 0) = JdrA [/li(r) + ;/l 2 (r)] 


(34) 


The values of qa for different as and different quark masses m are shown in Fig. [6] (for details 
of the numerical calculation see the Appendix). At = 1.23 (which provides a good £t to 
QCD phenomenology) we get qa = 0.63, which is close to the value obtained in Ref. |^. The 
axial charge exhibits only a minor dependence on the parameters of chiral symmetry breaking. 
So we can conclude that the discrepancy with phenomenological value is a qualitative feature 
and can not be cured by a simple tuning of the model. We note that this result is similar to 
the Skyrme model, where also has a small value 33(], which is weakly affected by the pion 
mass. The problem of small qa is a long standing problem in Skyrmion physics. There might be 
signihcant N^. corrections to this value, which are not taken into account here and may req uire 
a qualitative modihcation of the model. For a detailed discussion of this issue see e.g. |mI, l32| 
and references therein. 
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g. 



FIG. 6: The baryon axial charge for various us and various quark masses: from top to bottom 
mzm = {0.01,0.02,0.04,0.08} = 323MeV). 

VII. CONCLUSION 


In this paper we have studied the interplay between chiral symmetry breaking and the 
baryon in a holographic QCD setup. Although the numerical calculations were done for a 
specific “hard-wall” model, our qualitative treatment is fairly general and the results should 
hold for a generic model. We have proved in a clear-cut manner the conjecture advocated in 


Refs. 
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26| that there are essentially two regimes in baryon physics where the dependence of 


the baryon mass on the chiral condensate is very different. The numerical calculations show 
that the physical value of the baryon mass is in the transition region between the two regimes, 
which clearly demonstrates that the Ioffe’s formula for the baryon mass [n| derived in QCD 
sum rules is relevant. 

The main finding of our work is the phenomenon of repulsion of the baryon from the IR wall 
due to the interaction with the scalar field dual to the chiral condensate. It is very important to 
understand how this repulsion looks from the perspective of the Sakai-Sugimoto setup. At large 
values of the chiral condensate the baryon is moved away from the tip of the cigar. The only 
way to do this without violating P-symmetry is to split the soliton into two symmetric halves. 
One of them is moved to the Left brane and another to the Right one. For symmetry reasons, 
these two halves carry equal half-integer baryon charges and opposite half-integer Abelian axial 
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charges, producing in total the solution with unit baryon charge. 

In spirit the splitting of the holographic baryon is very similar to the splitting of an instanton 
on compact space into monopoles. This phenomenon takes place when the system is considered 
at finite temperature or finite chemical potential. However, the details can be very different 
and deserve a thorough study. In particular it is not completely clear how the splitting occurs 
in terms of the branes. The baryon is represented by the D4 brane wrapped around the internal 
sphere and cannot be naively split into half D4 branes. The process could be viewed upon a 
T-duality transformation around the compact coordinate on the cigar. Another delicate issue 
concerns the behavior at the nonvanishing temperature. Different scenarios of the splitting 
phenomenon could take place once one takes into account the temperature dependence of the 
chiral condensate. 

We should recall, that once one moves from the tip of the cigar, the separation between D8 
and D8 branes grows with the radius of the compact circle r. In our treatment we neglected 
this separation in favor of the opportunity of working with the 5D model with local fields, but 
in a general setup once the baryon is split and repelled from the tip of the cigar this separation 
may lead to interesting phenomena. For example, although the net Abelian axial charge of the 
solution is zero, the axial dipole moment is nonzero and proportional to the size of the r-circle. 
The 9- dependence of the Skyrmion due to the splitting can be considered in the SS model as 
well since the 6'-term corresponds to the holonomy of the one-form RR field along the compact 
coordinate on the cigar. One should note though, that the splitting effect is proportional to 
the size of the instanton and therefore is parametrically suppressed by the’t Hooft coupling in 
the original Sakai-Sugimoto model. On the other hand, the divergent tachyonic mode of the 
open string can provide a sufficiently large scalar boundary value, which will render the effect 
sizable. Given these subtleties, it is not clear in which particular setup one should study the 
phenomena mentioned above. 

The other interesting finding, that we made in the present study, is the backreaction of the 
baryon solution on the local density of the quark bilinear operator. We find that the chiral 
condensate gets suppressed inside the baryon and the chiral symmetry is partially restored. 
This idea was widely discussed in the literature but to our knowledge this phenomenon has 
not been directly observed in any models of QCD. The chiral condensate in the SS model 
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follows from tachyon condensation on the string, hence it would be interesting to investigate 
the impact of the Skyrmion on the tachyon action directly. This could also be related to the 
recent discussion of the spectrum of hadrons including baryons in terms of the details of the 


Dirac operator spectrum in QCD 
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Appendix: Accuracy of the numerical calculation 

We perform our numerical analysis using two different methods implemented in different 
codes in order to exclude any possible systematic errors (see text). The mass and topological 
charge distribution of the soliton is calculated on the grids 60 x 60, 120 x 120 and 240 x 240. 
One can observe a good convergence of the results (see Fig. [ 7 ]) and the estimated precision of 
the hnest grid is within 1%. 

It is numerically challenging to calculate the axial charge (l3Tj) . The mean value of the axial 
current is dehned by the asymptotic form of the solution near the boundary y = e —)■ 0. The 
equations of motion are singular on this boundary, hence the precision of the calculated solution 
should be higher than e. On top of that, the integrand in Eq. fl3T)) has several factors of r and 
in order to get hnite results one needs to ensure that the fall off of the solution at large r is 
calculated with a precision of order 1/r^. An additional difficulty is the fact that we impose 
nontrivial boundary conditions at r ^ oo. Because of that, the accuracy of the solution at 
r —)■ cxD is limited by the precision of the hnite difference approximation of the scalar held 
derivatives on the grid. At the end of the day we are forced to use very dense grids in the 2 ; 
direction in order to obtain convergent integrals for qa- The sample of the integrand in Eq. (134)) 
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FIG. 7; The accuracy of calculated baryon mass on different N x N grids with N = {60,120,240} 
(from top to bottom). is the extrapolated value for N ^ oo. 

for a = 2.6 and m = 0.01 calculated on the grids 60 x 120, 60 x 240, 60 x 480 and 60 x 1000 
is shown in Fig. [HI While for intermediate values of x, the results of all the grids coincide, the 
tail at a: = 1 (r —)■ oo) is resolved only on the hnest one. 


dgA 



FIG. 8: The integrand ([341) calculated for cr = 2.6 on different grids with Ny = {120,240,480,1000} 
(from top to bottom). The tail is decently resolved only on the finest grid. 
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